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We present exact analytical solutions to study the coherent interaction between a single photon and the 
mechanical motion of a membrane in quadratic optomechanics. We consider single-photon emission and 
scattering when the photon is initially inside the cavity and in the fields outside the cavity, respectively. 
Using our solutions, we calculate the single-photon emission and scattering spectra, and find relations 
between the spectral features and the system's inherent parameters, such as: the optomechanical coupling 
strength, the mechanical frequency, and the cavity-field decay rate. In particular, we clarify the conditions 
for the phonon sidebands to be visible. We also study the photon-phonon entanglement for the long-time 
emission and scattering states. The linear entropy is employed to characterize this entanglement by treating 
it as a bipartite one between a single mode of phonons and a single photon. 

The hybrid coherent coupling 1 between electromagnetic and mechanical degrees of freedom is at the heart of 
cavity optomechanics 2 4 . In general, optomechanical couplings can be classified into two categories: linear or 
quadratic couplings. Namely, the coupling term is proportional to either x or x 2 (x being the mechanical 
displacement). For a mechanical resonator, a linear coupling corresponds to a force acting upon the mechanical 
resonator and this leads to a displacement of its equilibrium position. However, for a quadratic coupling, it will 
change the resonant frequency of the mechanical resonator (in the new representation, rather than in the original 
representation). This is because the quadratic-coupling term can be integrated into the potential energy of the 
harmonic oscillator (changing the effective stiffness of a spring), and hence the frequency of the mechanical 
resonator is renormalized. This renormalized harmonic oscillator is related to the original one by a squeezing 
transformation. 

To better understand and exploit optomechanical couplings, it is highly desirable to realize these couplings in 
the single-photon strong-coupling regime, in which the couplings involving a single photon can produce 
observable effects on both mechanical and electromagnetic signals. Such a regime is important to test the 
fundamentals of quantum theory 5 " and to explore possible applications of optomechanical devices to future 
quantum technology 911 . In the past several years, much attention has been paid to the single-photon strong- 
coupling regime of linear coupling 12 25 . Considerable theoretical studies, such as photon statistics 1617,2123 and 
mechanical-state engineering 24 , have been carried out in this regime. Two theoretical proposals 26,27 have recently 
been suggested to reach this regime using superconducting circuits with Josephson junctions. In addition, 
experimental advances in linear optomechanics are being made towards the single-photon strong-coupling 
regime 1215 . However, for the quadratic coupling, though much attention has been paid to this area 28 4S , not 
much work has been devoted to the single-photon strong- coupling regime because the currently attainable 
coupling strength is weak. Recently, some methods have been proposed to increase the quadratic optomecha- 
nical coupling strength or to seek other possible realization of quadratic optomechanics. For example, an 
experiment 42 demonstrated that the quadratic coupling strength can be increased significantly using a fiber 
cavity with a smaller mode size, and a smaller and lighter membrane. A measurement-based method has also 
been proposed to obtain an effective quadratic optomechanics 36 . In addition, some other systems, such as 
trapped cold atoms or a dielectric nano- or microparticle, have been suggested to simulate an effective quadratic 
optomechanical coupling 44 . These works provide a possibility of studying the quantum nonlinearity in quadratic 
optomechanics. Motivated by these advances, it is of interest to study the quadratic optomechanical coupling in 
the single-photon strong-coupling regime. 

When a quadratic optomechanical cavity works in the single-photon strong-coupling regime, the frequency 
change of the mechanical resonator induced by a single photon will, in turn, significantly affect the cavity field, 
causing some observable features in the cavity photon spectrum. Thus, a natural question arises: how the 
spectrum may characterize the single-photon strong-coupling regime? In this paper, we answer this question 
by calculating analytically the spectrum of single-photon emission and scattering. In particular, we build a 
connection between the spectral features and the system's inherent parameters. We also clarify the condition 
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for observing the phonon sidebands in the spectra. It should be 
pointed out that below we assume that the mechanical resonator 
has been pre-cooled to a low-phonon-number regime. Moreover, 
in typical optomechanical systems, the optical decay rate y c is much 
larger than the mechanical decay rate )>m- Under these two condi- 
tions, we will only consider the optical dissipation and neglect the 
mechanical dissipation. This is because the emission and scattering 
processes will be completed in a time interval l/y c <^t<^l/y M . 
During this period the mechanical dissipation is negligible. In this 
issue, the emission and pulse-scattering cases are different from the 
continuous-wave driving case. When the system is driven by a con- 
tinuous-wave field, the mechanical dissipation should be included 
for the steady-state solution 46 . 

Accompanying the processes of single-photon absorbtion and 
emission, the total system experiences transitions involving phonon 
sidebands, and hence the frequency of the emitted photon will be 
related to the states of the phonon sidebands due to energy conser- 
vation. This relation leads to the generation of photon-phonon 
entanglement. Since the emitted photon exists in the continuous 
modes of the outside fields, this entangled state involves a single 
mode of phonons and a set of continuous modes of a single photon. 
In general, it is hard to characterize such type of entanglement. 
However, from the point of view of a single photon rather than 
photon modes, we could treat this entanglement as a bipartite one 
between a single photon and the phonon mode. For a pure initial 
state, the linear entropy can be employed to describe this 
entanglement. 

Results 

Quadratic optomechanical system. We consider a quadratic opto- 
mechanical system with a "membrane-in-the-middle" configuration 
[see Fig. 1(a)], where a thin dielectric membrane is placed inside a 
Fabry-Perot cavity. We model the moving membrane as a harmonic 
oscillator and focus on a single-mode field in the cavity. When the 
membrane is placed at a node (or antinode) of the intracavity 
standing wave, the cavity field will quadratically couple to the 
mechanical motion of the membrane. Let us denote the position 
and momentum operators of the membrane as x and p, then the 
Hamiltonian of the system is (with fi = 1) 

£> 2 1 

H ovc = m c (x)(aa+l/2) + — + -MQ 2 % 2 , (l) 
opc c\ j\ ' ' 2M 2 

where a ' and a are the creation and annihilation operators of the 
single-mode cavity field, respectively. In the quadratic coupling case, 
the cavity-field frequency depends on the mechanical motion by 

, , , 1 d 2 co c (x) 

co c {x) = oj c + qxr, with r\ = - — — - — | I= 0) where co c is the cavity- 
field frequency when the membrane is at rest. The parameters M and 
Q in Eq. (1) are the mass and frequency of the mechanical mode. By 
reorganizing the coupling term between the zero-point energy of the 
optical mode and the mechanical motion into the mechanical 
potential energy, we have 

p 2 1 1 

H opc = mca'a + + -Ma> 2 M x 2 + ria^ax 2 + - co c , (2) 



where ojm = yfl 2 + »;/M is the renormalized frequency of the 
membrane. We note that it is justified to work in the 
representation of the mechanical frequency oj m , because the 
coupling term rix 2 /2 always exists. By introducing the mechanical 
creation and annihilation operators V and b, by 
x = v /l/(2Mco M )(fo t + b) and p = i v / Moj M /2(& t - b) , the Hamil- 
tonian, up to a constant term (co c + oj m )/2, becomes 28 

H 0 p C = w c a' a + (OMb i; b+goa ; a(b' ! + £>) , (3) 



(a) Coupling mirror Membrane Perfect mirror 




Photon number 



Figure 1 | Quadratic-optomechanical system and the energy-level 
structure, (a) Schematic diagram of a quadratic optomechanical system 
with a "membrane-in-the-middle" configuration, (b) The diagram of the 
energy-level structure (unsealed) of the optomechanical system when the 
cavity is in a vacuum or contains a single photon. 



where g 0 = ri/(2Moj M ). The third term in Eq. (3) describes a 
quadratic optomechanical coupling with a strength g 0 between the 
cavity field and the membrane. We point out that since the 
reorganized frequency (Dm of the mechanical mode also depends 

on the coupling strength go by com = \J Q 2 + 2g 0 a>M, the ratio of 
the coupling strength g 0 over the mechanical frequency a> M is 
bounded by g 0 /oj M < 1/2. However, below we also consider 
parameters by extending the scope of the ratio beyond this bound. 
This is because the bound can be exceeded in some other quadratical 
optomechanical systems, such as trapped cold atoms. 

When there are s photons in the cavity, the last two terms in 
Hamiltonian (3) can be renormalized as a harmonic-oscillator 

Hamiltonian with the resonant frequency \J com ( ojm + 4sgo ) . To 
keep the stability of the membrane (i.e., the frequency should be a 
positive number), the strength g 0 should satisfy the condition (co M + 
4sg 0 ) > 0. The photon number operator a f a in Hamiltonian H opc is a 
conserved quantity, and hence for a given photon number s, the 
coupling actually takes a quadratic form sg 0 (b f + b) 2 , which can be 
diagonalized with the single-mode squeezing transformation. 
Denoting the harmonic-oscillator number states of the cavity field 
and the membrane as \s) a and \m) b (s, m = 0,1, 2, ...) respectively, 
then the eigensystem of the Hamiltonian H opc can be obtained as 

H 0 ^s)Jm(s)) b =(so Jc + SU + mco { ll)\s) a \m(s)) h , (4) 

where we introduce the s-photon coupled membrane's resonant fre- 
quency a>M and energy-level shift <5 <s) , 
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co$ = co My Jl+^, «M = i(a,W-« aM ). (5) 
The s-photon squeezed phonon number state in Eq. (4) is defined by 

\m(s)) b = S b (t,^)\m} b , (6) 



where Sbyi/^J =exp 
the squeezing factor 



V( fo2 - fot2 ) 



f/ s, = -ln 1 + 

4 V oj M 



is a squeezing operator with 



We note that the eigensystem of the Hamiltonian (3) has been 
derived in previous studies 31,45,46 . In the zero-photon case, we have 

I" 1 ("))(, = \ m )b' m< M =cu Af> and (5 <0) = 0. For following convenience, 
the energy-level structure of the system in the zero- and one-photon 
cases is shown in Fig. 1(b). 

To include the dissipation of the cavity field, we assume that the 
cavity photons can couple with the outside fields through the coup- 
ling mirror. Without loss of generality, we model the environment of 
the cavity field as a harmonic-oscillator bath. Then the Hamiltonian 
of the whole system including the optomechanical cavity and the 
environment can be written as 



H = H opc + j m k c[c k dk+^ f [ac[ + c k <$ ) dk, 



(8) 



where the annihilation operator c k describes the kth mode of the 
outside fields with resonant frequency co k . The coupling between 
the cavity field and the outside fields is described by the photon- 
hopping interaction with strength £. Since the decay rate y M of the 
mechanical resonator is much smaller than the decay rate y c of the 
cavity, then during the emission and wave-packet scattering time 
interval 1 / y t .<0<C 1/}'m> the damping of the membrane is negligible. 
In this work we take into account the dissipation of the cavity and 
neglect the mechanical dissipation. 

States in the single-photon subspace. In the rotating frame with 
respect to H 0 = w c a^a + co c I cj,c k dk, the Hamiltonian (8) becomes 



Hi = co M tfb+gocfia(tf + b) + \ A k c' k c k dk 

J o 

+ £/ (acl + c k a')dk, 



(9) 



where A k = w k — co c is the detuning of the kth mode photon from the 



cavity frequency. The total photon number, N = a~ a + I c k c k dk, in 

J o 

the whole system is a conserved quantity because of [N, Hj] = 0. 
Denoting \m(l)) b = \rh) b for conciseness, a general state in the 
single-photon subspace of the total system can be written as 



|p(f)>= ^A m (f)|l)» 



00 p GO 

+ E 

m = 0 



I B m , k (i)\0} a \m} b \l k )dk, 
J o 



(10) 



where |l)Jw) fc |0) denotes the state with the cavity in the single- 
photon state |l) a , the membrane in the single-photon squeezed 
number state \rh) b (hereafter we call it as squeezed number state 
for conciseness), and the outside fields in a vacuum |0). Also 
|0)al w )i|lfc) denotes the state with a vacuum cavity field |0) a , the 
membrane in the number state \m) b , and one photon in the kth 



mode of the outside fields \ l k ). The variables A m (t) and B mk (t) are 
probability amplitudes. 

We point out that these squeezed number states in Eq. (10) satisfy 
the completeness _ fl \ m) b h (m\—I b (I b is the identity operator in 
the Hilbert space of mode b) and orthogonality b{™\n) b = & m , n - 
Moreover, the overlap b {m\h) b = b(m\S b yl^j \ n )b between the 
squeezed number state \n) b and the harmonic-oscillator number 
state \ m) b is determined by the relation 

, — ; — r Floor [^j Floor jjj 



(7) b{m\S b {n (1) )\n 



(cosh^ 1 ))"^ 2 



E E 

f = 0 / = 0 



(-!)' 
/!/'! 



Q tanhi? 



(i) 



(11) 



-21)1 



cosh rf 



J m-2!,n-2h 



where r\ m = (1/4) ln(l + 4g 0 /oj M ) and the function Floor[x] gives the 
greatest integer less than or equal to x. In principle, the state of the 
whole system can be obtained by solving the Schrodinger equation 
under a given initial condition. Below we will consider single-photon 
emission and scattering. 

Single-photon emission. In the single-photon emission case, a single 
photon is initially inside the cavity, and the outside fields are in a 
vacuum. Without loss of generality, we assume that the initial state of 
the membrane is an arbitrary number state | n 0 ) b . Once the solution in 
this case is obtained, the solution for the general initial membrane 
state can be obtained accordingly by superposition. In this case, with 
the Laplace transform method, we obtain the long-time solution for 



these probability amplitudes as A n 



= 0 and 



R , , ^ ^jW) b (^\n)bb(^ 0 ) b e-'^ + ^ 

B na , m ,k{<X))= 2^ " "Tfj m — — . (12) 

„=o Aji-(5 l ' —ncoy +ma)M + iy c /2 

where we introduce the cavity-field decay rate y c = 2n(; 2 and add the 
subscript n 0 in A„ 0>m (t) andB„ 0imi A;(r) to mark the membrane's initial 
state \ n 0 ) b . 

In the long-time limit, the single photon completely leaks out of 
the cavity and hence the cavity is in a vacuum [A„ o m (co) =0]. The 
amplitude 5„ o m j t (co) exhibits a clear physical picture for the single- 
photon emission. Specifically, the initial state |l) B |«o)fc ca n be 
expanded as E^o c «o.«|l)«l")l>> with c n a ,n = b(n\n 0 ) b . For each 
component |l) fl |w) fc , the single-photon emission process induces 
the transition |l) a |»)(, _> |0) a |m) ft . The corresponding transition 
amplitude is proportional to the numerator in Eq. (12). Due to the 
quadratic terms of b and b 1 in S b (if i] ), \m) b and \n) b should have the 
same parity, i.e., being odd or even. Consequently, the phonon num- 
ber distribution in the long-time state of the membrane will have the 
same parity as its initial component j n 0 ) b . In addition, we can derive 
the resonant condition in this emission process from the energy-level 
structure in Fig. 1(b). For the transition |l) II |w) (j ->|0) IJ |m) fc , the fre- 
quency of the emitted photon is m k = (O c + d^ + *"-■• 11 



which is consistent with the resonance condition 



Ajt = <5 (1) + — mco M , 



(13) 



obtained from the pole of the denominator in Eq. (12). 

We know from Eqs. (10) and (12) that, corresponding to the initial 
state j l) a |«o)(,|0)> the long-time state of the whole system is 



koM)=E/ B n,,m,kM\0) a H b \h)dk. (14) 
m = 0 ■' 0 

Therefore, when the membrane is initially in a general density matrix 
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(0)= £ P&(°)l»>»»<»l. 



US) 



the long-time state of the whole system can be obtained by super- 
position as 



P (00)= J2 P^)\ Vm (co))( 9n (co)\. 



(16) 



Here the density matrix elements are p„ n [0) = b{m\p (0)|«) fe . 

To characterize the quadratic optomechanical coupling, a useful 
quantity is the single-photon emission spectrum, i.e., the probability 
distribution of the emitted photon. For the initial membrane state 
(15), the emission spectrum is defined by 

S(A k ) = Tr[\l k )(l k \p(co)] 

/,m,n — 0 

In Fig. 2, we plot the emission spectrum S(A k ) versus the photon 
frequency Aj-, for various values of g 0 and y c , when the membrane is 
initially in its ground state |0)(,. We see from Eq. (12) that both 

> 'Yc an d m M > Vc might be the resolved-sideband condition. 
For a positive g, then w M > y c could make sure that the two condi- 
tions are met, because of coj^' >o>m- We found that, in the case of 

0J ii' > Tc >cy M» the phonon-sideband evidence is negligible. So, in 
this paper, we consider co M > y c as the resolved-sideband condition. 
Figures 2(a-c) are plotted in the resolved-sideband regime co M > y c 
so that the phonon sideband peaks could be used to characterize the 
coupling strength go When g 0 < y c [Fig. 2(a)], the spectrum is 
approximately a Lorentzian function with width y c and center 
Afc = 8 m . In this case, there are no sideband peaks in the spectrum. 
However, the sideband peaks become visible when g 0 > y c . 
Physically, when the displacement of the membrane equals its 
zero-point fluctuation, the photon frequency shift induced by the 
quadratic optomechanical coupling is g 0 . To resolve this frequency 
shift from the Lorentzian spectrum of a free cavity, the condition 



(a)g> =0.1 



' 1 ' 1 - 

(b)g» / »M = 0 - 4 



A, 
2 




v M M' M 



k M 

Figure 2 | Single-photon emission spectrum S(A^) versus for various 
values of g 0 and y c . Panels (a-c) are plotted in the resolved-sideband 
regime (yJco M = 0.2), while (d) is plotted in the unresolved-sideband 
regime (yJoj M = 1.5). The membrane's initial state is I0)j,. 



g 0 > y c should be satisfied. Such a condition can also be understood 
by examining the height of these peaks in the spectrum. To resolve a 
peak in the spectrum, the peak height should be much higher than 
the tail of its neighboring Lorentzian. This requires go~S>)' c in the 
resolved-sideband regime. As an example, we analyze the special case 
ofgo/cojvf^l- m the resolved-sideband regime com/Vc^I an d under 
the initial state |0)i,, we expand S(A k ) up to second-order in g 0 lco M . 
Then, the height of the sideband peaklocated at A k = <5 (1) — 2oj m can 

be obtained as s(<5 (1) -2co M J ~ (y c /&noj 2 M ) (l + 8gg/y 2 c ). Since the 

main peak of the spectrum is approximately a Lorentzian function 



Si (A, 



271 



-<5 (1) ) 2 +^/4 
s(s {1) -2wmJ 



, then the requirement 



SdS® 



+ ^r» 1 

2(D M ) Yc 



(18) 



leads to the condition go^>}' c - 

We remark that the positions of these sideband peaks in Fig. 2 are 
determined by the resonance condition (13), and these sideband 

peaks are not periodic because of the difference between u)^ and 
<x>m- Also, for the initial state |0)f, of the membrane, the contributing 
m and n in Eq. (13) should be even numbers due to the parity 
requirement. In Figs. 2(a-c), the peak located at A k = <5 (1> is the main 
peak (corresponding to the transition | l) fl |6) (j — > |0) a |0) (j ). Hence, we 
can resolve the main peak from the Lorentzian spectrum for a free 
cavity when (5 (1) > y c , which requires g 0 > y c (l + yJoj M ). The peak 
located at A^ = S m — 2a> M corresponds to the transition 



|1) B O). -H0)j2) & . Moreover, the peak located at 



A k = S^ 



2cQm — 2ojm corresponds to the transition |l) a |2) (j ->|0) a |2) (; . 
Finally, Fig. 2(d) is plotted in the unresolved-sideband regime. We 
can see from Fig. 2(d) that, even though the system works in the 
single-photon strong-coupling regime, there are no sideband peaks 
in the spectrum. This fact indicates that the resolved-sideband 
regime and the single-photon strong-coupling condition are two 
combined necessary requirements for observing sideband peaks in 
the emission spectrum. 

To illustrate how the spectrum depends on the initial state of the 
membrane, we plot in Fig. 3 the spectrum S(A k ) versus A k when the 
membrane is initially in either the Fock states |0)f, and |l);„ the 
coherent state \ ji = 1)/,, or the thermal state pf(n = 1). For the initial 
state |0) fc , the main peak (with the location A k = <5 (1> ) in Fig. 3(a) is 
related to the transition | \) a jo}^— > |0) a |0) fe . The two peaks located at 



<5 (1) - 2co M and A k = 5 {1) +2w 



(i) 



■ M —2com correspond to the 
transitions |l>„|6) fc -*|0>j2)j and |l>«|5> fc -*|0>.|2> 6 , respectively. 

main peak (located at 



For the initial state |l)b, the 



A* = <5 
|1).|1>»- 



►|o>ji; 

(i) 

M ' 



-CO 



-<x>m) in Fig- 3(b) is related to the transition 
b . The other two peaks located at 

-3com and A k = <5^ + 3cu^' — 3ojm correspond to 



the transitions H>,|l)j-»|0>,|3>,, and |l) fl |3} i> ^p) B |3) fcJ respect- 
ively. Here, Figs. 3(a) and (b) only show even- and odd-parity side- 
band peaks, respectively. However, the coherent and thermal states 
contain both odd- and even-parity number states, and hence we can 
see both odd- and even-parity sideband peaks in Figs. 3(c,d). The 
positions of these sideband peaks are consistent with those in 
Figs. 3(a,b). 

In the long-time limit, though the single photon is completely 
leaked out of the cavity, its state is still entangled with the mechanical 
mode. This entanglement involves a single mode of phonons (the 
mechanical degree of freedom) and a set of modes of the photon 
because the single photon is distributed into the continuous fields 
outside the cavity. In general, it is difficult to clearly describe the 
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). 0) - W 

+2 co -2 &>,,)/ co, , 

M M' M 



(^ 1) +3cu'' ) -3ffl„)/ ( y„ 

r M M' M 



< 
Co 



0.5- 




k M 

Figure 3 | Single-photon emission spectrum S(A t ) versus \ k for various 
initial states of the membrane. The membrane is initially prepared in 
either the Fock states 10);, and I the coherent state I ft = 1);,, or the 
thermal state p t t(n= 1). Other parameters are yJ(o M = 0.2 and go/w M = 
0.6. 




Figure 4 | The linear entropy £jin the emission case. The linear entropy £; 
versus the optomechanical coupling strength g 0 , when the initial state of 
the membrane is 10);,, 1 1)^, and 12);,. Here yJ(a M = 0.2. 



corresponding to the membrane's initial states 1 0>^» | l)fc> an d |2)b, the 
linear entropy keeps increasing. We may roughly explain this phe- 
nomenon by analyzing the magnitude distribution of these factors 
b(n 0 \s') b and b (s\n 0 ) b . When n 0 changes from 0 to 2, the magnitude 
distribution of these transitions elements b {n 0 \s') b (for different 5') 
becomes increasingly smoother. This implies that the number of 
contributing coefficients increases, and hence the entanglement 
increases. 



structure of this entanglement. However, from the viewpoint of a 
single photon, we can characterize the entanglement as a bipartite 
one between a single photon and a single mode of phonons. In 
particular, we will consider a pure initial-state case so that the 
long-time state of the total system is also pure; then we can employ 
the linear entropy to quantity this bipartite entanglement. 

When the membrane is initially in the general state (15), the long- 
time state of the total system is given by Eq. (16). In terms of Eq. (12), 
the reduced density matrix of the membrane can be obtained as 



pW(oo) = % 



E pU°) ne 



-i(l—l')<OMt 



/,/',m,n,s,5' — 0 



where 



b{l\s} bb {~sH bb {n\~s'} bb {~s'\l'} b 
(1-1')oj m + 0' - s) com + iy c 



\l) bb (ll (19) 



(20) 



The linear entropy 49 of the density matrix (19) is 

1 2^ 



£; = 1-Tr 



E 



E pUm 



l,l' = 0 m,n,s,s' = 0 



(21) 



In Fig. 4, we plot E t versus g 0 f° r initial Fock states \n 0 ) b of the 
membrane, i.e., p h m (0) = <5 mi „ 0 <5„, no . Figure 4 shows that (as a general 
trend) Ej increases with increasing g 0 . However, there are some res- 
onance dips in the linear entropy when g 0 takes some special values. 
The locations of these dips can be determined from the poles of the 

denominator in Eq. (20), i.e., (I — 1')ojm + (s' — s)co^ =0, here the 
values of (I — V) and (s' — s) should be even numbers because of the 
parity requirement in the transitions. For example, there is a dip atg 0 / 
a>M = 0.75, which corresponds to | (Z — Z')/(s' — s)\ = 2. In addition, 



Single-photon scattering. In the single-photon scattering case, the 
single photon is initially in a Lorentzian wave packet 

/zr 1 , , 

t dk in the outside fields, where A 0 and £ are the 

a A k -A 0 + k 

detuning center and spectral width of the photon. In this case, with 
the Laplace transform method, the long-time solution (i.e., t2>l/y c , 
1 /e) of these probability amplitudes is obtained as A„ 0im ( 00 ) = 0 and 



1 

71 A,t — A 0 + ie 
1 



(22) 



nA k - [A 0 + (n 0 - m)w M ] + it 

x J2 i ycb(m\n) bb (h\n 0 ) b \ 

„ =0 A k - (5 (1) - nfflj^ + mco M + iy c /2J 

Here the subscript n 0 in these amplitudes is used to mark the initial 
state of the membrane. In the long-time limit, the single photon will 
completely leak out of the cavity, and hence we have A„ o m (co) = 0. It 
can be seen from B„ o m ^(co) that there are two physical processes in 
the single-photon scattering, (i) The single-photon direct- reflection 
process: the incident photon is directly reflected by the mirror, 
without entering the cavity. This process is described by the first 
term of B„ 0jm ^(co). (ii) The photon-membrane interacting process: 
the single photon enters the cavity to couple with the moving 
membrane and eventually leaks out of the cavity via the cavity 
decay channel. This process is described by the second term (i.e., 
the second and third line) in Eq. (22). In this process, the system 
experiences the transitions |0) I! j«o)(, — > !l)al")fc — > l")a! m )(i' These 
transitions are governed by the two resonance conditions 



A 0 = <5« + 



nm^ -n 0 co M , 



-5® 



(1) 

+ nar M ' — mco M , 



(23a) 



(23b) 
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Figure 5 | Single-photon scattering spectrum S(A t ) versus \ k for various 
values of g 0 and y a Figures (a-c) are plotted in the resolved-sideband 
regime (yJco M = 0.2), while (d) is plotted in the unresolved-sideband 
regime (y c /co M = 1.5). The membrane's initial state is 10);,. Other 
parameter are A 0 = d m and e/oj M = 1.2. 

which can be derived from either the energy level structure in 
Fig. 1(b) or the poles of the probability amplitude (22). 
Interestingly, the second line in Eq. (22) is a Lorentzian wave 
packet with spectral width e and center = A 0 + (n 0 — m)w M . 
In comparison to the initial Lorentzian wave packet, the shift of the 
wave packet center is equal to the energy variance of the membrane. 
Moreover, the third line in Eq. (22) has a similar form as Eq. (12) for 
the single-photon emission process. 

The single-photon scattering spectrum can be calculated in terms 
of Eqs. (17) and (22). We see from Eq. (22) that either the second line 
or the third line could cause phonon sidebands, and the conditions 
for resolving these sidebands due to the two lines are com > £ and oj m 
> y c , respectively. Here, y c is the system's inherent parameter while e 
is an externally controllable parameter. In the following, we first 
consider the case of £ > u>m > y c so that the observed sideband peaks 
are caused purely by the system inherent effect. In Fig. 5, we plot the 
spectrum S(Ajt) versus the photon frequency A^ for various values of 
g 0 and y c . When g 0 < y c , there are no peaks in the spectrum 
[Fig. 5(a)]. The phonon sideband effect can be observed in the scat- 
tering spectrum when g 0 > Yc and <J) M > y c [Figs. 5(b,c)]. Owing to 
the interference between the direct reflection process and the 
photon-membrane interacting process, there exist both peaks and 
dips in the spectrum. In Figs. 5(b,c), the dips represent the transition 
l-Oal^);, - > |0)al")ii» while the peaks correspond to the transition 
|l) a |5), — >|0)j2) 6 . In addition, in the unresolved-sideband regime 
(y c > u>m), there are no peaks even in the single-photon strong- 
coupling regime [Fig. 5(d)]. 

We now consider the near-monochromatic case (£<Cy c ). In 
Fig. 6(a), we plot the scattering spectrum in the case of y c > co M 
and £<^u>m- This figure exhibits phonon sideband peaks, and hence 
indicates that <±>m > £ also provides the condition for observing the 
phonon sideband peaks due to the second line in Eq. (22). We point 
out that this provides a way to characterize the coupling strength g 0 
from the scattering spectrum in the case of y c > co M . Another benefit 
in the near-monochromatic case is that we can conveniently control 
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Figure 6 | Single-photon scattering spectrum S(A t ) versus \ k for various 
A 0 and initial states of the membrane. The panel (a) is plotted in the 
unresolved-sideband regime (go/oj M = 2 and }'J(o M = 1.5), while panels 
(b-d) are plotted in the resolved-sideband regime (go/oj M = 0.8 and }'Joj m 
= 0.2). The frequency center A 0 of the incident photon is 6^ +2oj^', and 
<5 a) in other panels. The membrane's initial state is the ground state 1 0);, in 
panels (a-c) and the coherent state a = 1)^ in (d). The parameter 
e/i>m = 0.02. 

the exciting transition by choosing the frequency of the incident 
photon. In Figs. 6(b,c), we choose the frequency of the incident 

photon as A 0 = d (1) and A 0 = <5' 1 ' +2oj^ to resonantly excite the 
system from |0)„j0) b to 11)^6)^ and |l) a |2),, respectively. In the 
emission process, the membrane will experience the transitions from 
15}^ and [2)^ to |w)j, (n = 0, 2, 4, ...). Therefore, the maximal 
frequency sideband peaks should be located at A^ = <5 (1) and 



2m^, respectively. In addition, the period of these peaks is 
2co M . Similar to the emission case, the scattering spectrum also 
depends on the initial state of the membrane. In Fig. 6(d), we plot 
the scattering spectrum when the membrane's initial state is coherent 
state I a = 1)^. Though the initial coherent state contains both even- 
and odd-parity states, the spectrum only exhibits similar peaks as 
those in Fig. 6(b). This is because the incident photon (with A 0 = <5 (1) ) 
only resonantly excites the membrane from |0);,to 6) fe ; other transi- 
tions from \n) b to \n') b (n, n' = 1,2, 3, with the same parity) are 
significantly suppressed due to the large detuning. A further photon 
emission process induces the transitions from 1 6} , to |n)j, (n = 0, 2, 4, 
...). 

Similar to the emission case, the scattered photon is completely 
emitted out of the cavity in the long-time limit; and the state of the 
photon is entangled with the mechanical membrane. For the mem- 
brane's initial state |n 0 )fc, the linear entropy of the long-time reduced 
density matrix of the membrane was calculated exactly. This is not 
shown here because the analytical solution is long. We can examine 
how the linear entropy depends on the system parameters: the incid- 
ent photon frequency A 0 and the optomechanical coupling strength 
go- In the near-monochromatic limit c<^y c , we plot, in Fig. 7(a), the 
linear entropy £; as a function of A 0 when the quadratic optomecha- 
nical coupling go takes various values. We can see from Fig. 7(a) that 
there are resonant peaks in the entropy. For a given g 0 , the locations 
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Figure 7 | The linear entropy Ej in the scattering case, (a) The linear 
entropy £; versus A 0 for various values of g 0 . The initial state of the 
membrane is 10};,. (b) The linear entropy £; versus g 0 when the membrane 
is initially in states 1 0);„ I and 1 2);,. Here the frequencies of the incident 
photon are A„ = <5 (1 >, A 0 = d^ +co^ — c»m> and A 0 = (5 (1) +2co^' — 2ca M , 



respectively. Other parameters are yJa>M = 0.2 and t/ra M 



"At 
= 0.02. 



of these peaks are A 0 = <5 (1) and A 0 = <5''' + 2ar^ , which are deter- 
mined by the resonance conditions in the dominant transitions in the 
photon injection process: |0)JO) i ,^|l) a |5) 6 and |0)j0) fc ^|l) a |2) fc . 

We also investigate the dependence of the linear entropy £; on the 
coupling strength go when the incident photon is in resonance with 
the transitions. In Fig. 7(b), we plot the entropy Ej versus g 0 when the 
initial state of the membrane is |0)f„ |l);„ and |2)(,. Here the single 
photon is resonantly injected into the cavity. Corresponding to the 
initial states |0);„ |l){,, and |2)f„ the driving frequencies are A 0 = 

A 0 = + (0"m — com, and A 0 = <5* 1 ' +2oj^' — 2co M , respectively. 
These drivings determine the dominant photon-injection transitions: 
|0>.|0>»-»|l> a |0) i . |0)Jl),-|l) 8 |l) i , > and |0) a |2)^|l) 8 |2) ii . We 
can see that, in the resonant scattering case, the linear entropy 
increases when increasing g 0 . 

Discussion 

Although the currently- available quadratic couplings are too weak to 
reach the single-photon strong-coupling regime, advances have 
recently been made in the enhancement of this coupling strength. 
In quadratic optomechanics, the coupling strength is g Q = t]x^ { , 
where x zpf is the zero-point fluctuation of the mechanical membrane, 

1 S 2 0) c {x) , 

2 dx 2 

frequency. Recently, the value of 17 has been increased significantly 
from about 30 MHz/nm 2 (in Ref. [33]) to 20 GHz/nm 2 (in Ref. [42]) 
using a fiber cavity with a smaller mode size. For a x zp f —41 fm, 
suggested in Ref. 28, the coupling strength is g 0 ~ 2n X 5.35 Hz. If 



and i-j -- 



| x = 0 , with w c (x) being the x-dependent cavity 



the cavity decay rate y c ~ MHz, the coupling strength needs to be 
further increased by five orders of magnitude to reach the single- 
photon strong-coupling regime. From g 0 = fJXJpf , we can see that the 
go can be increased by obtaining a larger r\ or x zp f. This requires the 
improvement of experimental conditions. In addition, the model 
under investigation is a general quadratic optomechanical 
Hamiltonian 36 . It can also be realized in various physical systems 
such as ultracold atoms and superconducting circuits. In the linear 
optomechanical coupling case, the ultracold atom system has been 
demonstrated to approach the single-photon strong-coupling 
regime 13,14 , and the superconducting circuit system has been esti- 
mated to be in this regime 26,27 . Moreover, other methods have been 
recently explored to achieve effective strong quadratic optomecha- 
nics. For example, a measurement-based method has been proposed 
to perform this mission. Based on these achievements, it might be 
possible to pursue the single-photon strong- coupling regime in 
quadratic optomechanics. 

In the above discussions, we did not include the mechanical dis- 
sipation. Now we give a rough estimate for the influence of the 
mechanical dissipation on our results. In our considerations, all 
the results are determined by the probability amplitudes given in 
Eqs. (12) and (22). To evaluate the influence of the mechanical dis- 
sipation, we introduce an imaginary decay factor into the resonant 
frequency of the mechanical mode, i.e., approximately replacing cu M 
with co M — iy M (2n t h + 1)/2, where h t i, is the thermal phonon occu- 
pation. In this way we can estimate the effect of the thermal mech- 
anical dissipation. If y M (2n t i 1 + l)<Cy c /2, then the mechanical 
dissipation is negligible. This is because, in the low-temperature 
regime, the imaginary part in the denominators of Eqs. (12) and 
(22) can still be approximated by y c /2, and the exponential factor 
e -7A((2«,/, + i)t is amiost 1 during the time scale r ~ lly c . 

We have analytically studied the single-photon emission and scat- 
tering in a quadratically-coupled optomechanical system. By treating 
the optomechanical cavity and its environment as a whole system, we 
have obtained the emission and scattering solutions using the 
Laplace transform method. Based on our solutions, we have calcu- 
lated the single-photon emission and scattering spectra, and found 
relations between the spectral features and the system's inherent 
parameters. In particular, we have clarified the condition under 
which phonon sideband peaks can be observed in the photon spectra. 
In the resolved-sideband regime cu M > }' c > the phonon sidebands are 
visible when g 0 > y a while the condition for resolving the photon- 
state energy-level shift <5 (1> is g 0 > y c (l + yJto M ). We have also 
investigated the creation of photon-phonon entanglement in the 
emission and scattering processes. This entanglement was created 
due to the energy requirement in the photon absorption and emis- 
sion processes, and it was treated as a bipartite one from the view- 
point of a single photon rather than photon modes. We have 
considered the pure state case so that the linear entropy can be used 
to characterize this entanglement between the phonon mode and a 
single photon. 

Finally, we want to mention several possible applications of the 
single-photon emission and scattering processes, (i) Since the emis- 
sion and scattering spectra are related to the eigen-energy of the 
system, we can read out the system parameters, such as coupling 
strength and mechanical frequency, from the spectra, (ii) We can 
infer the initial state information of the mechanical resonator based 
on the single-photon spectra, similar to the linear optomechanics 
case 50 . As a special example, the initial temperature of the mechanical 
membrane can also be determined from the spectra, (iii) The entan- 
glement generated in the emission and scattering processes involves 
phonons and one outgoing photon. So this entanglement could be 
used to implement quantum information processing tasks, where 
phonons and photons are used as information memory and carriers, 
respectively. 
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Methods 

Solving the equations of motion for probability amplitudes. According to Eqs. (9), 

d 

(10), and the Schrodinger equation t — |(p(f)) =Hj\<p(t)}, we obtain the equations of 
motion for the probability amplitudes 

A m (0 = -*(<5 (1) +ma)<^ (24a) 



K,k(t) = -i(A k + ma M )B m *(t)-i£^ b {m\n) b A 1t (t), (24b) 

ji = 0 

where the single-photon coupled membrane's frequency and energy-level shift are 
given by oJj$ = m M y/l + 4g 0 /oj m and 3^ = [&m ~ 0J M^j j 2- In addition, the 
coefficients b(m\n) b and \,{m\h) b can be calculated using Eq. (11). 

The equations of motion (24) for the probability amplitudes may be solved with the 
Laplace transform method under a given initial condition. In the single-photon 
emission case, a single photon is initially inside the cavity, and the outside fields are in 
a vacuum. For the mechanical mode, we first assume that its initial state is an arbitrary 
number state |hq)j. Based on the solution in this case, the solution for the general 
initial membrane state can be obtained accordingly by superposition. For the initial 
state |l)Jwo)fe|0)> tne corresponding initial condition for Eq. (24) is 
A m (0) —b{fh\no) b and B mt k(0) = 0. In the single-photon scattering case, the single 
photon is initially in a Lorentzian wave packet in the outside fields, and the cavity is in 
a vacuum \0) a . We also assume that the membrane is initially in the number state 
|«o)&» anQ then the initial condition for Eq. (24) becomes A m (0) — 0 and 

B m fc(0) — \ — t S m n , where A 0 and c are the detuning center and spectral 

V 7i A k — A Q + ic 

width of the photon, respectively. Based on these initial conditions and the equations 
of motion, we may obtain the solution of these probability amplitudes using the 
Laplace transform method. The long-time solutions corresponding to the single- 
photon emission and scattering have been given in Eqs. (12) and (22), respectively. 
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